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ABSTRACT 

Soft phase lags, in which X-ray pulses in lower energy bands arrive later than pulses in higher 
energy bands, have been observed in nearly all accretion-powered millisecond pulsars, but their origin 
remains an open question. In a study of the 2.5 ms accretion-powered pulsar SAX J1808.4— 3658, 
we report that the magnitude of these lags is strongly dependent on the accretion rate. During the 
brightest stage of the outbursts from this source, the lags increase in magnitude as the accretion rate 
drops; when the outbursts enter their dimmer flaring-tail stage, the relationship reverses. We evaluate 
this complex dependence in the context of two theoretical models for the lags, one relying on the 
scattering of photons by the accretion disk and the other invoking a two-component model for the 
photon emission. In both cases, the turnover suggests that we are observing the source transitioning 
into the "propeller" accretion regime. 

Subject headings: binaries: general — stars: individual (SAX J1808.4— 3658) — stars: neutron - 
stars: rotation — X-rays: binaries — X-rays: stars 



1. INTRODUCTION 

Accretion-powered millisecond pulsars (AMPs) pro- 
vide a unique laboratory for stud ying the process of disk 
accre tion onto magnetic stars (jPsaltis fc Chakrabarty! 
1 19991 and refrences therein). The ~10 8 G stellar mag- 
netic field truncates the accretion disk only a few neu- 
tron star radii above the surface, channeling the infalling 
matter toward the magnetic poles. The resulting X-ray- 
emitting hot spots and accretion shocks are modulated 
by the stellar rotation to produce the observed pulsa- 
tions. Such compact accretion geometry poses a partic- 
ularly interesting challenge. We want to understand the 
details of the interaction between the magnetic field and 
the disk, and the effects that result from the proximity of 
different emitting regions. The accretion-powered 2.5 ms 
pulsar SAX J1808.4— 3658 is a particularly good system 
for this type of study: in addition to having multiple out- 
bursts, it spans orders of magnitude in accretion rate. In 
this paper we report our study of the energy-dependence 
of pulsations from this source and show that it gives new 
insight into these problems. 

SAX J1808.4-3658 wa s the first detected AMP 
(|Wiinands fe van der K lis 1998), and it remains the 
best studied. Persistent pulsations have been detected 
throughout the four ~1 month outbursts from this source 
that have been observed by the Rossi X-ray Timing 
Explorer (RXTE) during 1998-2005. The timing of 
these pulsations establishes a 2.01 h r binary orbital pe- 
riod (|Chakrabarty fc Morgan! 119981 ) that is increasing 
on a timescale of P or b/-forb = (66 ± 4) x 10 6 yr, pos- 
sibly due to the ablati on of its low-mass c ompanion 



during X-ray quiescence (Hartman et al.l 2008. hereafter 



H08; iDi Salvo et "all I2008D . The presence of X-ray pul- 
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sations during the peak accretion rate and the long- 
term spin-down of this source provide a tight constraint 
on its magnetic field: B = (0.4-1.5) x 10 s G (H08; 
iPsaltis fc Chakrabartvlll999f) . A distance of 3.4-3.6 kpc 
is estimated from thermonuc lear X-ray burst properties 
([Galloway fc Cumminjl2006f ). 

The outbursts from SAX J1808.4-3658 had similar 
light curves, which we divide for convenience into four 
stages (see Fig. 2 of H08 for a schematic): a rise 
from quiescence lasting «5 days; a short-lived peak 
at a 2-25 keV flux of (1.9-2.6) x 10~ 9 erg cm" 2 s -1 ; 
a 10-15 day slow decay in luminosity that lasts until 
the source reaches «0.7 x 10~ 9 erg cm~ 2 s^ 1 ; and fi- 
nally a rapid fall in luminosity followed by weeks of 
low-luminosity flares with a ~5 d periodicity, during 
which the observed f lux can vary by a factor of 1000 
(|Wiinands et al.ll200lD . This sudden drop has been inter- 
pret edas_atransjti£nJnto a "propeller" accretion regime 
(e.g.. lGilfanov et al . 1998) i n which accretion is only par- 
tially inhibited (| Spruit fc Taaml 119931 : iRappaport et all 
|2004[) . although it may also be caused by a change in 
the viscosit y of the accretion disk as it transitions in to 
quiescence ( Gil fanov et al]|l99& iGierlihski et"a^[2002T ). 

Shortly after t he discoy e ry o f pulsations in 
SAX J1808.4-3658, iCui et all (fl998h noted the pres- 
ence of soft phase lags, the tendency of the pulses 
at lower ("softer") energies to arrive later than the 
corresponding pulses at higher ("harder") energies. The 
magnitude of the lags increased approximately linearly 
between 2-10 keV and saturated at 0.08 rotational 
cycles (200 /xs) above 10 keV. Similar soft lags have 
been detected in six other AMPs: IG R J00291+5934 
(iGallowav et al.ll2005l ). X TE J0929-314 dGalloway et all 
I2002D. XTE J175 1-305 (IGierlihski fc Poutanenl l2005h. 
XTE J1807-294 (iKirsch et al.l l2004D. XTE J1814-338 
([Watts fc StTohmaver1l200i r~and HETE J1900.1-2455 
(jGallowav et al.ll2007h . 

A number of explanations for the AMP soft lags 
have been proposed. Models that include o nly Doppler 
boosting, due to the stellar rotation ([Fordl 120001 : 
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IWeinberg et al.l |2001[) or the orbiting m aterial in the 
accretion disk (jSazonov fc Sunvaevl 120011), have short - 
comings that are summarized bv lGierlihski eTaT](|2002T) . 
The two most plausible explanations invoke a two- 
component emission model with differing angular depen- 
dence and spectra (jPoutanen &: Gierlinskil 12003: here- 
after PG03) or the scattering delay of hot photons 
from the a ccretion shock off the coo ler disk or stel- 
lar surface (jFalanga k. TitarchuMl2007l hereafter FT07). 
In this paper we present an energy-resolved analysis 
of the persistent pulsations for all the outbursts of 
SAX J1808.4-3658 detected by RXTE. We will discuss 
the implications of our results for these models of the soft 
lags. 

2. DATA ANALYSIS 

For our study of the energy dependence of the per- 
sistent pulsations, we analyzed the same RXTE obser- 
vations of SAX J1808.4-3658 as in H08 (refer to Ta- 
ble 1 therein). The E_125US_64M_0_1S dat a mode of the 
RXT E Proportional Counter Array (PCA: IJahoda et all 
1996) was used for nearly all the observations. It pro- 
vides 122 fxs time bins and 64 energy channels divided 
evenly over the detector response, a good combination for 
discerning the energy dependence of the persistent pul- 
sations. The energy band covered by each channel varied 
over the 7 yr of observation spanned by t his study due 
to ad justments to the PCA anode voltages (|Jahoda et al.l 
2006); we accounted for this effect throughout our mea- 
surements. 

For our timing analysis, we shifted the photon arrival 
times to the solar system barycenter, applied the RXTE 
fine clock correction, and filtered out data during Earth 
occultations and intervals of unstable pointing. We also 
removed any data within 5 min of thermonuclear X-ray 
bursts. To measure the phases and fractional amplitudes 
of the persistent pulsations for a given energy band E, we 
folded 512 s intervals of data using the phase timing solu- 
tion derived in H08. Dividing the resulting pulse profiles 
into n phase bins with photon counts XE,j, j = 1 ■ • -n, 
the amplitude (AE,k) and pulse phase offset (<j>E,k) f° r 
the kth harmonic is 

n 

A E ke ^k4, E , k =2J2 ^je 2 " jl/ " • (1) 
i=i 

We measure the phase offsets of the energy bands rel- 
ative to the softest band, Eq: AtpE.k = 4>E,k — 4>E ,k- 
Thus negative values of A<f>E,k indicate that pulsations 
of the fcth harmonic within band E arrived earlier than 
the pulsations within the softest band. The uncertainties 
of these phase offsets due to Poisson noise are 



2-KkA E ,k 

where Ne is the total number of photons in the energy 
band (H08). For our analysis, we rejected phase mea- 
surements with <7E,fe > 0.04 cycles (i.e., 0.1 ms). 
The measured fractional rms amplitudes are 

rE ' k = V2(N E -B E y (3) 

5 Throughout this paper, we number the harmonics such that the 
fcth harmonic is k times the frequency of the 401 Hz fundamental. 
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Fig. 1. — Average phase offsets A4>E,k (relative to the 2.1—3.3 
keV band) and fractional amplitudes during the slow-decay stages 
of the 2002 outburst (black; MJD 52564.1-52575.0) and the 1998 
outburst (gre y; MJD 5 0919.0-50932.8, the same pointed observa- 
tions used by Cui et al. 1998). Negative offsets indicate soft lags. 

where Be is the approximate number of background 
events within our energy band and time interval, esti- 
mated using the FTOOL pcabackest. 6 Uncertainties on 
the fractional amplit udes are com puted using the method 
described bv iGrothl (| 1975( 1 and IVaughan et all (| 1994( 1. 
which correctly accounts for the addition of noise to the 
complex amplitude of the signal. 

We estimated the 2.5-25 keV fluxes of the persistent 
emission by fitting the spectra extracted from Standard2 
mode data (128 energy channels, 16 s time resolution) 
with an absorbed blackbody plus power law model, in- 
cluding an iron lin e if there were signifi cant residuals 
around 6.4 keV; see lGallowav et alJ (pool . §2.1, for full 
details. 

3. RESULTS 

To assess the energy dependence of the pulse phase and 
fractional amplitudes, we calculated the average phase 
offsets and fractional amplitudes during the 11 d slow- 
decay stage of the 2002 outburst. Figure Q] shows the 
results, along with measurements from the 1998 outburst 
for comparison. We chose this time range because it was 
the longest and most data-rich interval during which the 
pulse profile of SAX J1808.4— 3658 was observed to be 
reasonably stable. The pulse profiles during the slow- 
decay stage of the 2005 outburst, although similar to the 
2002 profiles, were somewhat less stable (see Fig. 3 of 
H08). 

The soft phase lag of the pulsations is immediately ap- 
parent for the fundamental. Between 2 keV and 10 keV, 
the pulses in harder bands arrive progressively sooner, 
with a slope of -0.014 ± 0.002 cycles keV -1 . Between 

6 http : //heasarc . gsfc.nasa.gov/docs/xte/recipes/pcabackest . html 




Fig. 2. — Phase offsets and fractional amplitudes as a function of energy for the four SAX J1808.4— 3658 outbursts observed by RXTE. 
The top panels show the background-subtracted light curves for each outburst. The strips along the tops of the graphs indicate the 
times of observations; stars indicate the times of the seven observed thermonuclear bursts. The dotted line indicates the critical flux of 
0.7 X 10 — 9 erg cm -2 s — 1 , at which the source transitions into the flaring-tail stage of the outburst. The second and third panels from the top 
show the phase offsets of the medium energy band (green; roughly 5—10 keV) and hard energy band (blue; £5:10-20 keV) pulsations relative 
to the soft band (red; £b2-5 keV). Negative phase offsets indicate that higher-energy pulses lead the soft pulses. The bottom two panels 
show the background-subtracted fractional amplitudes for these three energy bands. Each point in the phase and fractional amplitude plots 
gives a day-long average. (They are not evenly spaced at 1 d intervals due to the uneven spacing of the RXTE observations.) Data within 
5 min of thermonuclear bursts were excluded from these averages. 



10 keV and 20 keV, this trend saturates at a lead of 
0.08 cycles (200 fis). There is impressive agreement be- 
tween the 1998 and 2002 outbursts, which also had nearly- 
identical light curves; the phase lags during the dim- 
mer 2005 outburst are similar in shape but plateau at 
a smaller magnitude. The phase offsets of the second 
harmonic are less pronounced, disappearing almost en- 
tirely or becoming hard lags when the source is bright, 
as is the case for the data in Figure [TJ 

In contrast, the energy dependence of the fractional 
amplitudes varies considerably between outbursts. The 
fractional amplitude of the fundamental during the slow- 
decay stage of the 1998 outburst (5.2-7.5% rms) was con- 
siderably higher than during the 2002 outburst (3.8-5.0% 
rms), as can be seen in Figure [1] the amplitude of the 
2005 outburst (2.0-2.6% rms) was lower still. The mor- 
phology of the energy dependence differs as well. The 
1998 outburst exhibits a large fractional amplitude peak 
at 3.5 keV and falling amplitudes up to 20 keV, similar 
to wh at is seen in XTE J1814-338 (|Watts fc Strohmaverl 
120061 ). In contrast, during 2002 and 2005 the fractional 
amplitude vs. energy curves are much flatter, with small 
but significant peaks at 8 keV. Given the strong similar- 
ities between the light curves of the 1998 and 2002 out- 
bursts, it is curious that the energy dependence of their 
fractional amplitudes is markedly different. All three 



outbursts show weak evidence of the fundamental's frac- 
tional amplitude increasing again above «20 keV. The 
fractional amplitude of the second harmonic generally in- 
creases with energy for all the outbursts, indicating less 
sinusoidal pulse profiles in harder energy bands. XMM 
observations of the 2008 outburst of this source reveal 
that the fractional amplitude of both harmonics contin- 
ues to decrease below 2 keV (jPatruno et al.l l2009'). 

Figure [2] illustrates how the energy-dependent phases 
and fractional amplitudes change over the course of the 
four observed outbursts. The soft lags in the fundamen- 
tal are apparent throughout all four outbursts; for the 
second harmonic, soft lags are apparent everywhere ex- 
cept during the outburst peaks. However, the degree 
of these lags is not constant: they are much more pro- 
nounce d during the tails of the out bursts, when the flux 
is low. IWatts fc Strohmaverl (|2006f ) noted a similar pat- 
tern in XTE J1814— 338, although that source did not 
trace as wide a range of fluxes, so the relationship was 
less clear. The energy dependence of the fractional am- 
plitudes does not exhibit any clear correlation with the 
source flux or any other measured parameters, apart from 
the general trend of fractional amplitudes increasing in 
the outburst tails (see §3.4 of H08). 

To further understand the relationship between the 
soft lags and the flux, we measured the phase vs. en- 
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Fig. 3. — Correlation between the slope of the phase lags and 
the flux in SAX J1808.4— 3658. Each slope is determined by a 
linear fit to the phase offset vs. energy measurements over 2-10 
keV, with each point representing aid average. Large negative 
slopes indicate large soft lags. The data from all four outbursts 
are shown: 1998 as diamonds, 2000 as triangles, 2002 as squares, 
and 2005 as circles. The dashed lines indicate the best linear fits 
above and below the 0.68 X 10" 9 erg cm" 2 s" 1 turnover. These 
fits are based on the 2002 and 2005 data only, as the 1998 and 2000 
outbursts exhibit greater scatter. 



ergy slopes over 2-10 keV. We chose this range because 
it has the most significant energy dependence for the 
fundamental and, to a lesser extent, for the second har- 
monic. Figure [3] plots the resulting slopes as a func- 
tion of flux. For the fundamental, the magnitude of 
the soft lags is clearly anticorrelated with the 2-25 keV 
flux when the source is bright: as the source dims from 
its peak at 2.7 x 10 -9 erg cm" 2 s" 1 down to roughly 
0.7 x 10~ 9 erg cm" 2 s" 1 , the magnitude of the soft lags 
grows. This range of fluxes covers the peak and slow- 
decay stages of the outburst. At lower fluxes, the relation 
reverses, and the magnitude of the lags becomes corre- 
lated with flux. The flux at which this turnover occurs 
is significant: once it falls to 0.7 x 10~ 9 erg cm" 2 s" 1 , 
the source dims rapidly and enters the flaring-tail stage 
of the outbursts. 

The trend for the second harmonic is somewhat less 
clear. Measuring a slope through the 2-10 keV fractional 
amplitude points does not work as well since the energy 
dependence is less linear (top panel of Fig. QJ, so there 
is more scatter. Even so, the data strongly suggest a 
similar relationship at high fluxes, particularly when the 
relatively scattered 1998 data are excluded. 

The slopes of these trends differ significantly for the 
two harmonics. In the low-flux regime, the slopes are 



— 16 ± 2 for the fundamental and —2 ± 2 for the second 
harmonic; at high fluxes, they are 32 ± 4 for the funda- 
mental and 17 ± 3 for the harmonic. (The units in all 
cases are 10 -3 cycles keV -1 [log flux] -1 .) These slopes 
were fit using the 2002 and 2005 data only, as the two 
earlier outbursts had appreciably more scatter and less 
data. Note that these units reflect rotational cycles of 
the star, so a mechanism that induces the same phase 
lag in all photons of a given energy would result in the 
slopes being equal for the two harmonics. For the high- 
flux slopes, we can conclude with 8 a significance that 
this is not the case. 

We also looked for correlations between the lags and 
the spectral parameters from our blackbody plus power 
law fits. The only parameter for which a correlation is 
apparent is the normalization of the hard power law com- 
ponent. This relation is not surprising, because the hard 
normalization and the total X-ray flux are linearly pro- 
portional: during the 1998-2005 outbursts, log(TVhard) = 
(1.02 ± 0.03) log(/ a; ), with a reduced x 2 near unity. The 
uncertainties of -/Vhard are much greater than the flux un- 
certainties, particularly when the source is dim, and the 
large scatter obscures whether a turnover is present in 
the lag vs. Ah ar d relation. 



4. DISCUSSION 

The energy dependence of the pulsations from 
SAX J1808. 4-3658 exhibit strong variability both be- 
tween and during the four outbursts observed by RXTE. 
The overall fractional amplitudes and the shape of the 
amplitude vs. energy curves differ substantially between 
outbursts, with no obvious correlation with flux; the only 
clear trend is that the pulsations generally become less 
sinusoidal with energy in the observed 2-30 keV range. 
In contrast, the soft lags of the pulses have a consistent 
energy dependence, increasing in magnitude in the 2- 
10 keV band and saturating at 10-30 keV, with a strong 
dependence on flux. At the highest fluxes the lags are 
near zero, but they increase steadily as the flux falls. 
Below a 2-25 keV flux of w0.7 x 10" 9 erg cm" 2 s" 1 , this 
trend reverses and lags get smaller again. This particular 
flux also marks the point in the outburst where the lu- 
minosity decays rapidly and the source enters the flaring 
tail. 

Soft lags are a common feature of AMP pulsations, 
and two plausible models have emerged to explain their 
existence. PG03 suggested a two-component model in- 
volving soft emission from a surface hot spot and a 
hard Comptonized component from the shocked layer 
in the accretion funnel. Each component has a dif- 
ferent beaming pattern, which is affected in a differ- 
ent way by the Doppler effect. This results in the 
soft pulse lagging the hard pulse. While this model 
works well for SAX J1808.4- 3658 and XTE J1751-305 
(|Gierlinski fc Poutanenl [20051 ) . it does not predict the 
reduction in soft l ags at very h i gh en ergies seen in 
IGR J00291+5934 (jFalanga et alJ 120051 ) . As a result 
FT07 proposed an alternative model in which the disk 
plays a prominent role. In this model, the main pulsed el- 
ement is the hard component from the accretion column: 
this radiation is then downscattered by cold plasma in 
the inner disk (or on the neutron star surface) to gener- 
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ate a soft lagging component. 7 These soft lags get larger 
as the density of the scattering plasma falls. 

Let us now explore the factors that might cause the 
flux dependence of the soft lags, assuming that flux di- 
rectly traces accretion rate. We then need to ask how the 
properties of the system are expected to change as the 
accretion rate falls. For a star with a magnetic field that 
is strong enough to channel the accretion flow, the reduc- 
tion in accretion rate should cause an increase in inner 
disk radius (e.g., iGhosh et al.lll977[) . This will have two 
relevant effects. First, the density of scattering plasma 
near the star will fall: within the model of FT07 this 
would cause the soft lags to increase. In addition the 
relative angular velocity of the inner disk and star will 
reduce. This will affect the an gle at which the accre - 
tion funnel impacts the surface (|Romanova et al.ll2004j ). 
changing the column density above any surface hot spot. 
This type of change might work within the model of 
PG03 to affect the soft lags. 

If these changes continue steadily, however, it is hard 
to see what could cause the turnover seen in Figure [3) 
One intriguing possibility that would give a very nat- 
ural explanation for such a turnover is that the sys- 
tem reaches the point where the magnetospheric radius 
r m exceeds the corotation radius r c . This is often re- 
ferred to as the "propeller" regime, since it was origi- 
nally thought that matter would then be ejected from 
the sy stem, with the star going from spin-up to spin- 
down (jlllarionov fc Sunvaevlll975T ). In fact, the transi- 
tion is probably not sharp, and r m need s to be substan- 
tially larger than r c for ejection to occur (| Spruit fc Taaml 
1993). At intermediate accretion rates, continued ac- 
cretion is possible even while the neutron star starts t o 
spin down (ISpruit fc Taamlll993l:lRappaport et~al"1l2004f ) . 
Simulations bv lLong et all ( 20051 ) indicate that the tran- 
sition between spin up and spin down should occur when 
r m /r c 0.7 for AMPs. The observed turnover happens 
at an accretion rate of M = 1.8 x 10~ 10 Mq yr _1 , giving 
Tnil^c = 0.6 for the magnetic field B = 1.5 X 10 8 G from 
H08. 8 This close agreement between the simulations 
and the M at which the flux rapidly drops and the soft 
lags turnover strongly suggests that SAX J1808.4— 3658 
passes through this transition as the accretion rate falls. 
In this eventuality there are two effects that could cause 

7 Ibragimov & Poutancn (2009) recently argued that photoelec- 
tric absorption in the accretion disk would prevent the large num- 
ber of scattering events required by the FT07 model. However, 
stand ard disk models give a tem perature of ~300 eV in the inner 
disk (Shakura & Svunvacv 1973); the actual temperature will be 
higher due to the irradiation. At these temperatures and the ex- 
pected electron densities (n e ~ 3x 10 19 cm -3 ; see FT07 and Fig. [4] 
here), all elements through oxygen will easily be fully ionized by the 
hard tail of the thermal distribution. The only significant source 
of photoelectric absorption is Fe, but for this n e only Fe xxv and 
higher remain at T > 300 e V, with the transitio n to full ioniza- 
tion occuring at 400-600 eV l|Rubiano et al.1120041) . Its He-like and 
H-like forms have K edges at 9 keV, so they will not absorb the 
softer photons responsible for the lags. Since the FT07 model is 
not ruled out a priori, this discussion attempts to address both 
models without preference. 

8 Accretion rates were estimated in the usual manner, by equat- 
ing the bolometric and infall luminosities: iitd? f x Cho\ = GMM/R, 
using d = 3.5 kpc and the bo lometric correction cj, i = 2.12 from 
IGallowav &~ Cumming ( 2006]). The Alfven radius was estimated 
as r m = kA^GM )- 1 / 7 M~ 2 / 7 (BR 3 ) 4 ? 7 , with k A = 0.5 as per 
I Lone et all 1 120051) . M = 1.4 Mq and R = 10 km were used in all 
cases. 
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Fig. 4. — Electron density n e of the inner accretion disk, as 
fit using the model of FT07, versus the mass accretion rate for 
the 2002 and 2005 outbursts of SAX J1808. 4-3658. The vertical 
dotted line marks the M at which the magnitude of the soft lags 
is greatest and below which the outburst is in its flaring-tail stage; 
a transition from normal accretion to the "propeller" regime can 
explain this change in the observed behavior. The dashed lines 
give the best fits of n e oc M 1 in each region. 

the turnover in soft lags. 

The first is due to the fact that the relative angular ve- 
locity of the inner disk and the star changes sign. The ac- 
cretion funnel will go from being dragged forwards to be- 
ing dragged backwards (assuming that it settles rapidly; 
iRomanova et al.l I2004T ) . The column depth above any 
surface hot spot should therefore go through a maximum 
near the zero torque point, causing a turnover in the scat- 
tering properties of the accretion column. If the accre- 
tion column plays a major role as discussed in PG03, this 
could cause the observed turnover. Another explanation 
within the PG03 model is the suppression of the lagged 
soft component at high fluxes due to the increased optical 
depth of the accretion shock when M is large. However, 
this mechanism would require the normalization of the 
soft component to decrease during the outburst peaks, 
contr ary to what is seen (cf. Fig. 2 of iGierlihski et all 
[2002T) . 

A second possibility stems from the fact that the prop- 
erties of the disk must also change at this point. For ac- 
cretion to continue after zero-torque point, the inner edge 
stops retreating, and the d ensity at the inner edge must 
rise (|Spruit fc Taaml 119931 ). In the scenario of FT07, 
this rise in inner disk density would cause the magni- 
tude of the soft lags to decrease again. In the absence of 
Compton upscattering, which would introduce hard lags 
not observed in SAX J1808.4— 3658, their model predicts 
that the scattering from the disk will modify the arrival 
times as At(E) « (m e c/aTn e )(E^ 1 — E~ r ), where ctt is 
the Thompson cross-section and n e is the electron den- 
sity of the scattering region. This model provides rea- 
sonable fits to the observed phase lags; Figure [4] shows 
the resulting n e . 

The fitted electron densities probe the relationship 
between the inner disk and the accretion rate. n e ~ 
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10 19 ' 5 cm -3 gives a mean free path of £ = (<7T«-e) 1 ~ 
0.5 km and time £/c ~ 2 /is, requiring ^100 scatter- 
ing events to produce the observed soft lags. This I 
rules out the thin neutron star atmosphere as a scat- 
tering medium. On the other hand, simulations by 
iRomanova et a l. (2004) suggest an inner disk thickness 
of order the neutron star radius, and (R/£) 2 ~ 400 is 
more than adequate to accomodate the number of needed 
scattering events. From the fits in Figured! we find that 
M 10±01 during the main outburst; during the pre- 
sumed propeller stage, n e oc ilif - °" 4±0 * 1 . An interesting 
consequence of the former relationship is that at higher 
accretion rates before the turnover, when r m /r c < 0.7, 
the flow rate (the volume of material per unit time cross- 
ing the inner edge of the disk) must be constant: since 
M = flmpner^Vr, where fl is the solid angle of the inner 
edge and v r is the radial velocity of material crossing it, 
the flow rate £lr^v r is fixed. 

If the phase lags in SAX J1808.4— 3658 are entirely due 
to Compton downscattering, as suggested by FT07, then 
this model must explain why the magnitudes of the soft 
lags for the second harmonic are roughly half the mag- 
nitudes for the fundamental over a wide range of M, 
as can be seen in Fig ure [3] Following the formalism of 
ISunvaev fc Titarchukl (|1980f ). we model the timing ef- 
fects of Comptonization with a distribution P( At) of the 
photon escape times from the scattering medium. For 
the fcth harmonic of the spin frequency v a , we convolve 
the harmonic component x k (t) = Ak exp(2mkv s t) with 
P(At) to find the new amplitude A' k and the time lag At 
of the scattered pulses: 



A I 2Kikv s {t—At) _ 



A k e 



2-nik.Vs t' 



P(t' - t)d£ (4) 



In general, the choice of P{At) accounts for the geome- 
try of the emitting and scattering regions. However, for 
the large number (~100) of scatterings previously noted, 



the geometry becomes unimportant and a photon diffu- 
sion disWbjrdon_of_P(A^ ^r" 1 exp(Ai/r) is appropri- 
ate (jSunvaev fc Titarchukj [T980). Applying this form to 
equation ((4]), the time lag of the k th harmonic is 



At k = (2irkv s T) 



_1 tan 1 27rfc^,r . 



(5) 



The diffusion timescale cannot be longer than the lags 
themselves: r < Ai max = 200 fis. This sets a lower 
limit of At%/Ati > 0.85 due to scattering alone, contrary 
to the observation that Abz/Abx w 0.5. Thus Compton 
scattering cannot by itself fully explain the observed lags: 
a change in the energy dependence of the initial pulses 
is still necessary. 

Clearly, detailed calculations are required to check 
whether the two models really would predict the out- 
comes that we suggest, but they are at least in principle 
plausible. One would also need to test whether the mod- 
els would predict any reversal in other pulsation prop- 
erties at propeller onset, given that no such changes are 
observed. 

The interesting behavior of the soft lags in this source 
suggests that the lags of the other accretion-powered mil- 
lisecond pulsars should be revisited to check for flux de- 
pendence. The only source for which this has been done 
is XTE J1814— 338: an incr ease in soft lags was observe d 
in the tail of the outburst (|Watts &; Strohmaveri r2006'). 
but the rapid drop off and shorter tail of this source 
made it very difficult to study the pulsations at the low- 
est fluxes. The other AMPs may be more promising as 
candidates for further study. 
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